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ABSTRACT 

In order to study the inner parts of the circumstellar material around optically faint, infrared 
bright objects, we present the first medium-resolution spectropolarimetric data taken in the 
near-infrared. In this paper we discuss Pa/3 line data of GL 490, a well-known embedded 
massive young stellar object, and of MWC 349A and MWC 342, two optically faint stars that 
are proposed to be in the pre-main sequence phase of evolution. As a check on the method, 
the classical Be star £ Tau, known to display line polarization changes at optical wavelengths, 
was observed as well. Three of our targets show a "line effect" across Pa/3. For £ Tau and 
MWC 349A this line effect is due to depolarisation by a circumstellar electron-scattering 
disk. In both cases, the position angle of the polarisation is consistent with that of the larger 
scale disks imaged at other wavelengths, validating infrared spectropolarimetry as a means 
to detect flattening on small scales. The tentative detection of a rotation in the polarization 
position angle at Pa/3 in the embedded massive young stellar object GL 490 suggests the 
presence of a small scale rotating accretion disk with an inner hole - similar to those recently 
discovered at optical wavelengths in Herbig Ae and T Tauri stars. 

Key words: polarisation - scattering - stars: circumstellar matter - early type - evolution - 
pre-main sequence 



1 INTRODUCTION 

Many open questions regarding the formation and evolution of stars 
require the ability to probe the inner parts of their circumstellar 
material, where accretion processes occur and stellar outflows find 
their origin. A powerful tool to address such issues is spectral line 
polarimetry. The method can utilise the fact that free electrons in 
ionised material scatter photons and thereby polarise photospheric 
radiation. In the case of spherical symmetry, all polarisation vec- 
tors cancel, and if the envelope is unresolved a net zero polarisa- 
tion is observed. On the other hand, if the circumstellar geometry 
is aspherical, such as a disk, a measurable net polarisation can be 
observed. This aspect of the method exploits the fact that ionised 
gas also emits recombination lines, which, by virtue of their lo- 
cation within the extended ionised gas undergo less scatterings by 
electrons, and will thus be less polarised. A "line-effect" is then 
observed. Electron scattering typically results in polarisation of or- 
der 1% (see Poeckert & Marlborough 1976) and occurs at scales 
of order several stellar radii (e.g. Cassinelli, Nordsieck & Murison, 
1987). The method can be used for objects that display emission 
lines, and immediately provides the answer to the question whether 
the circumstellar (electron-scattering) material is spherically sym- 
metric when the source itself cannot be resolved in imaging (e.g. 
Oudmaijer & Drew 1999, Harries et al. 2002, Schulte-Ladbeck et 
al. 1993). A further advantage of spectropolarimetry is that it pro- 
vides additional constraints to straightforward spectroscopy, which 



is particularly relevant for follow-up modelling (Drew et al. 2004, 
Harries 2000, Vink, Harries & Drew 2005a). 

Given the debate in the massive star formation community 
about the respective roles of isolated disk accretion versus com- 
petitive accretion in dense cluster environments or even mergers of 
lower mass protostars, it is crucial to add discriminating observa- 
tional constraints to this discussion. In particular, it has been ar- 
gued that the absence of accretion disks in massive young stellar 
objects (YSOs) is reason to question the relevance of the disk ac- 
cretion scenario (e.g. Wolfire & Cassinelli 1987; see also Norberg 
& Maeder 2000). Although a positive detection of a disk around a 
YSO should by no means be seen as proof for disk accretion, a con- 
tinuing absence of reliable disk signatures may eventually be taken 
as evidence against it. 

We have obtained spectropolarimetric data of the intermediate 
mass pre-main sequence Herbig Ae/Be stars using the Ha emis- 
sion line (Oudmaijer & Drew, 1999, Vink et al. 2002). More than 
half of the Herbig Ae/Be stars show the line-effect, while the po- 
sition angles of the circumstellar material are consistent with the, 
larger scale, disks that in some cases have been imaged by other 
techniques (see the compilation by Vink et al. 2005b). The high 
proportion of line effect detections in these objects suggests that all 
Herbig Ae/Be stars are surrounded by disks on small scales. This 
result provides indirect support to the notion that Herbig Ae/Be 
stars have grown via accretion. 

However, more massive stars with even stronger radiation 
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Table 1. Targets 



Name 


V 


./ 


tcxp (s) 


Pcont (%) 


©cont (°) 


GL 490 


14 


10.5 


4800 


15.59± 0.05 


130.1 ±0.3 


MWC 349A 


13 


7.0 


800 


3.60± 0.04 


176.6 ±0.3 


CTau 


3 


3.2 


240 


1.10± 0.005 


35.4 ±0.1 


MWC 342 


10.6 


~ 6 


320 


0.40± 0.03 


152.5 ±2.3 



Continuum polarisations are measured in line-free regions in the spectra. 
The error bars are internal errors, and small due to the large number of 
pixels that are averaged. The external errors are estimated to be 0.15% and 
2°. 

pressure have remained elusive. These most massive young stel- 
lar objects stay embedded in their parental cloud until after settling 
on the zero-age main sequence, which prevents the use of optical 
line polarimetry. It is therefore necessary to go to (near-)infrared 
wavelengths to study such objects. This has never been done so 
far, presumably because common-user polarisation optics at near- 
infrared optimised telescopes are rare. To do such a study, we first 
need to consider which line to observe, which is no trivial task. 
One has to take into account the intrinsic strength of the line, the 
large extinction towards these objects and the contribution of excess 
emission due to thermally re-radiating dust. At the longer A'-band 
wavelengths the sources suffer less from extinction, but the dust 
excess emission may set an upper limit to the observable line de- 
polarisations. This is because the magnitude of the "line-effect" is 
determined by the strength of the intrinsic continuum polarisation 
against which the change across the unpolarized line is contrasted. 
Any excess emission due to an extended envelope will effectively 
dilute the stellar continuum polarisation and consequently this re- 
quires more sensitive data. As these issues can not be answered 
easily pending dedicated observations, we decided to conduct a pi- 
lot study to observe optically faint objects and, for the first time, 
apply this method in the near-infrared. We chose to use the Pa/3 
line in the J-band, one of the strongest hydrogen recombination 
and at a wavelength where the observable continuum radiation is 
more likely to be dominated by the star. 

We observed four objects, MWC 349A, MWC 342 and GL 
490, and as a check £ Tau. The first two objects are optically faint 
and are proposed to be young pre-main sequence stars while GL 
490 is a well-known embedded massive young stellar object. £ Tau 
is a classical Be star, known to display a line effect at optical wave- 
lengths. It was observed as well to test and prove the concept. 



2 OBSERVATIONS & DATA REDUCTION 

The observations were taken during the night of the 26 th Septem- 
ber 1999 using the UKIRT telescope on Mauna Kea, Hawaii. The 
weather was reasonable, with a seeing of order 1.3 arcsec. The 
CGS4 spectrograph was used in conjunction with the 150 lines/mm 
grating and the long camera. The detector was a 256x256 InSb ar- 
ray which has a pixel size projected on the sky of 0.6 arcsec. During 
the observations, the slitwidth was 4 pixels (2.4 arcsec), and the slit 
was oriented East- West. A krypton arc-lamp was used for wave- 
length calibration. This set-up resulted in a wavelength range of 
1.256- 1.31 /im, with pixels of ~ 2. 1 x 10~ 4 ^im. The final spectra 
have a spectral resolution of 170 km s~ . The objects do not appear 
to be extended in our data. 



For the linear polarisation observations, the IRPOL2 instru- 
ment was employed. This polarimetry module consists of a MgF2 
prism and a half-waveplate. The QU spectra were recorded at var- 
ious positions of the waveplate, while the objects were nodded up 
and down the slit to facilitate sky-subtraction. The observing se- 
quence consisted of consecutive observations of source and sky (by 
nodding over the slit) for every position of the waveplate. Polarised 
and un-polarised standards were observed throughout the night. A 
log of the observations is provided in TableQ 

The data were reduced in IRAF, and included flatfielding, sky 
and bias subtraction, bad pixel masking and extraction of the spec- 
tra. The resulting spectra were then imported in the Starlink 1 Time 
Series and Polarimetry TSP software, in which the Stokes parame- 
ters were determined. Further manipulations were done in the Star- 
link POLMAP package for the polarisation spectra and IRAF for the 
intensity spectra. 

From the two unpolarised standards, HD 10476 and 
HD 202573, we determined the instrumental polarisation to be 
0.165% at an angle of 122.2°. This value was subtracted from the 
reduced data. The zero-point of the polarisation angle was deter- 
mined assuming the polarisation standard stars HD 23512 and HD 
43384 to have polarisation angles of 29.6° and 170.7° respectively. 
These values correspond to the R band values listed by Hsu & 
Breger (1982). The corrections were within 1° from each other. The 
polarisation measured in the standard stars at 1.28 fim is slightly 
(0.1%) larger than an extrapolation from their optical data using 
a Serkowski law would suggest. This is a well-documented effect, 
the polarisation in the near-infrared drops of shallower than in the 
optical (e.g. the review by Whittet, 1996). 

The errors in polarisation per pixel range from ~ 0.05 — 0.1% 
for the brighter stars to 0.25% for MWC 342 and MWC 349A and 
0.40% for GL 490. For the brightest objects, the error in the polar- 
isation is observed to increase with increasing wavelength. This is 
most likely due to low level fringing occurring in the optics. The 
final, external, accuracies are estimated to be of order 0.15% and 
2° respectively. 



3 RESULTS 

Here we discuss the results individually for the objects. We start 
by describing the results of the classical Be star £ Tau, and then 
continue with the massive young stars. In Table 1 the measured 
polarisations of the objects are presented. 

3.1 C Tau 

Because £ Tau has a long history of (spectro)polarimetric obser- 
vations, it is an ideal object to test the diagnostic power of near- 
infrared spectropolarimetry. Also, £ Tau is known to have a large- 
scale disk found through interferometry, consistent with the in- 
trinsic line polarisation at optical wavelengths (Quirrenbach et al. 
1997). Our data are presented in Fig.Q where the Stokes I (in- 
tensity) vector is plotted in the bottom panel, and the polarisation 
percentage and polarisation angle (PA) are displayed in the middle 
and upper panels respectively. The polarisation across the Pa/3 line 
shows a marked drop with respect to the continuum. 

The continuum polarisation of 1.1% is slightly lower and the 
PA is very close to the optical values of Quirrenbach et al. (1997, 
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Figure 1. The left hand plot shows the polarisation data of £ Tau represented as a function of wavelength in /im, the bottom panel shows the (continuum- 
normalised) intensity spectrum, while the middle and upper panels show the polarisation and position angle respectively. The latter two are rebinned to a 
corresponding accuracy in the polarisation of 0.05% . The right hand plot shows the Stokes QU vectors with the same binning applied. The excursion of the 
depolarisation in the QU plane implies a best fitting intrinsic PA of 32° . 



their mean values are about 1.5% and 31°). The small change in 
polarization from the optical to the NIR is not expected from a 
Serkowski law, where the dust polarisation drops off steeply with 
wavelength, but is consistent with the flatter dependence from elec- 
tron scattering. We can measure the intrinsic PA from the excur- 
sion across the line profile observed in the QU diagram; B = | x 
atan(AU/AQ) yielding 32±1°. As scattering in the optically thin 
case gives a PA perpendicular to the disk, this compares very well 
with the disk's position angle of— 58 ± 4° measured interferomet- 
rically by Quirrenbach et al. (1997). 

In summary, the near-infrared spectropolarimetry is a comple- 
mentary technique to optical spectropolarimetry and interferometry 
allowing us to move to longer wavelength ranges where we can ob- 
tain detailed observational constraints on types of object that would 
otherwise remain obscured. 



3.2 MWC 349A 

This well-studied object is mostly considered to be a massive pre- 
main sequence star (e.g. Danchi, Tuthill & Monnier 2001; Meyer, 
Nordsieck & Hoffman 2002) although it is sometimes proposed to 
be an evolved massive star (strong proponents are Hofmann et al. 
2002). It had not been previously observed by us at Ha as the object 
is too faint in the optical for medium resolution spectropolarimetry. 
Our data are plotted in Fig.|2| The Pa/? line is strongly in emission, 
whilst the He I 1.27 fim line, blueward of Pa/3, is clearly present 
as well. The observed continuum polarization of 3.5%, at a PA of 
~ 180°, is consistent with optical values (Yudin 1996). There is a 
clear signature of a line effect at Pa/?, particularly seen as a drop in 
polarization percentage across the line (middle panel). 

Representing the data in QU space, such as shown in Fig. [2] 
allows one to obtain the intrinsic PA of the polarisation indepen- 
dent of any foreground polarization, which only adds a constant QU 
vector to the intrinsic polarization. There is a strong effect over the 
line, with a magnitude of ~ 1.5% at a polarization angle of 21±1° 
as determined from a weighted least squares fit through the QU 
data points, as above. A similar, less strong effect at the same PA 



is visible for the helium line. There is one earlier report of a 'line- 
effect' in MWC 349A: Meyer et al. (2002) obtained low resolu- 
tion (7.5-10 A) optical spectropolarimetry, and reported strong line 
depolarisations. Earlier, broad-band optical polarimetry accompa- 
nied by narrow-band Ha polarimetry was obtained by Zickgraf & 
Schulte-Ladbeck (1989), but they did not report the line-effect, pos- 
sibly because it was smeared out over the large wavelength range 
covered by the Ha filter. 

We also note that although the polarization angle changes by 
about 10° over the emission line, this does not necessarily have 
to be an intrinsic property of the polarization. One measures the 
vector sum of both intrinsic and interstellar polarization. Due to 
this vectorial nature, a simple depolarisation can change into a more 
complicated profile with the addition of ISP (see e.g. Oudmaijer et 
al. 1998 for a marked change in the case of HD 87643). The linear 
excursion observed in QU space strongly suggests that we deal with 
a straightforward depolarisation. 

The obvious question that now arises is where does the scat- 
tering arise, in the well-known ionised North-South bipolar out- 
flow visible in the radio (White & Becker 1985) or the East- West 
circumstellar disk as is evident from imaging of hydrogen recombi- 
nation line masers (e.g. Planesas, Martin-Pintado & Serabyn 1992; 
Rodriguez & Bastian 1994)? As mentioned in the previous sub- 
section, in the optically thin case we will measure a polarization 
angle perpendicular to the disk plane, while modelling of Be star 
disks implies polarization percentages not much more than 2% 
(Cassinelli et al. 1987; Waters & Marlborough 1992), with the po- 
larization decreasing with inclination angle. 

The intrinsic position angle of the polarization is 21 ± 1°, 
which is almost perpendicular to the position angle of 107±7°, 
measured from the H30a line maser (Planesas et al. 1992) and the 
100±3°measured from interferometric images tracing warm dust 
(Danchi et al. 2001). This is consistent with the scattering occurring 
in a structure at the same PA. As the larger scale disk is observed 
to be close to edge-on, the depolarization of about 1.5% lends extra 
support to the idea that the spectropolarimetry, and associated line 
effects, are due to a circumstellar ionized disk. 
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Figure 2. As previous figure. The polarization data of MWC 349A rebinned to 025%. The dashed line is an enhanced version of the intensity spectrum 
intended to reveal the relatively strong, twice the continuum, He I A 1.27842/im line. The polarization spectrum shows an effect across both the Pa/3 and He I 
emission lines. The depolarization is smaller across the He I line because of the relatively larger contribution of the continuum to the observed polarization. 



Since MWC 349A is such a well-observed object, we now 
consider how our Pa/3 data fit in with existing ideas about the ion- 
ized gas around this object. Spectropolarimetry should be particu- 
larly helpful as it, in principle, probes the ionized gas close to the 
star. 



Hamann & Simon (1986, 1988) proposed a model of the cir- 
cumstellar material, based on velocity resolved optical and near- 
infrared spectroscopy of the star. They found that the line-width, 
and peak separation of resolved hydrogen recombination lines de- 
creased as function of distance to the star. The Helium recombi- 
nation lines showed the largest widths, and were suggested to be 
originating at the inner edge of a Keplerian rotating circumstellar 
disk at a distance of ~ 2 AU from the star. The hydrogen recom- 
bination lines, especially optically thick lines like Ha, were placed 
much further from the star, consistent with the observed extent of 
the hydrogen masers at 60 AU and more (e.g. Planesas et al. 1992). 
The authors explicitly introduced a low density region between the 
star and the disk/wind. This is mainly because of the observed line- 
widths: if MWC 349 is indeed a massive hot star, the line widths 
of order 100 km s~ x are too small to be identified with a hot stel- 
lar wind whose observed velocities are generally larger. Instead, 
Hamann & Simon locate the ionised gas in a photo-evaporating 
disk (see also Hollenbach et al, 1994), further away from the star, 
where the local escape velocity is low. This would present a coher- 
ent picture, where a comparatively low density fully ionized inner 
disk is responsible for the observed electron-scattering induced po- 
larization. 



We wish to mention a possible alternative that can explain the 
low outflow velocity and still be consistent with a circumstellar 
disk. Drew, Proga & Stone (1998) present a hydrodynamical model 
of radiation driven disk winds for YSOs that also predict narrow 
Hi line emission in the equatorial plane, and higher velocity, more 
tenuous emission in the polar regions. Such a model could also be 
responsible for the observed line-profiles in MWC 349A. 



3.3 MWC 342 

This object, often classed a young Herbig Be star, has a rich emis- 
sion line spectrum with Ha displaying a P Cygni profile, indicat- 
ing current outflow, while the presence of many forbidden emis- 
sion lines class the object a B[e] star (Jaschek & Andrillat 1999). 
The most in-depth study to date of MWC 342 is probably that of 
Miroshnichenko & Corporon (1999) who suggest the object is a 
luminous evolved star based on its extinction distance. The liter- 
ature does not report any observations that specifically study the 
geometry of the circumstellar material. Clues that at least some 
asymmetries are present come from the variable optical polariza- 
tion reported by Bergner et al. (1990). The continuum polarization 
around Pa/3 of this object is 0.4% This makes it the least polarized 
of the objects in our sample and corresponds to a phase of low opti- 
cal polarization (Bergner et al 1990, Miroshnichenko, private com- 
munication). The data are plotted in Fig.|5]but no line effect can 
be inferred with confidence. Since the polarization across the line 
measured per pixel is only at the 2<r level - as opposed to the value 
derived in Table Q which is derived from many pixels - it is hard 
to measure an effect at a significant level, and easy to overinterpret 
such data. Rebinning the data would yield severely undersampled 
data and potentially lead to misleading results. The uncertainty in 
polarization of the individual pixels is 0.25%, so the upper limit 
to the line-depolarisation is at least 0.75%, as an entire resolution 
element needs to be covered. 



3.4 GL 490 

AFGL 490 is a heavily embedded massive young star with a lumi- 
nosity of order 1O 3 L0 (e.g. Schreyer et al. 2002). It is subject to an 
optical extinction > 20 (Bunn et al. 1995), and is surrounded by a 
thick dusty disk at a PA of around 120-145 °. The disk structure has 
been observed using a variety of techniques: near-infrared imaging 
was done by Minchin et al. (1991); high resolution molecular CS 
observations were taken by Schreyer et al. (2002), whilst Alvarez 
et al. (2004) detect a, possibly transient, structure in near-infrared 
speckle interferometry. In addition, a flattened radio structure was 



Pa(3 spectropolarimetry 5 




found by Campbell, Persson & McGregor (1986). The central ob- 
ject powers a wide bi-polar CO outflow roughly perpendicular to 
the disk (e.g. Mitchell et al. 1995, Schreyer et al. 2002). 

The optical polarization was studied by Haas, Leinert & 
Lenzen (1992). They found that the bulk of the optical polariza- 
tion of ~18% can be explained by a combination of foreground 
polarization (surrounding stars have optical polarizations of up to 
13%) and polarization by circumstellar dust. As electron scattering 
typically would contribute around 1-2%, it is indeed most likely 
the dust scattering, rather than electron scattering that contributes 
most to the circumstellar polarization. Because the dust can not ex- 
ist close to the star owing to its comparatively low condensation 
temperature, the Pa/3 line forming region and the central source 
continuum will be seen as a point sources by the dust and be equally 
polarized. Although the observed polarization is a combination of 
interstellar dust and circumstellar dust and electron-scattering, the 
"intrinsic" electron scattered radiation is revealed in QU space, as 
the interstellar and circumstellar dust polarization will add respec- 
tive constant QU polarization vectors, whereas the intrinsic electron 
scattering polarization can be dependent on the velocities. 

Our polarimetric spectrum of GL 490 is presented in Fig. [4] 
The emission line itself is resolved as it is slightly wider than the 
instrumental profile, thereby allowing the study of changes in the 
polarization characteristics across the line itself. We measure a con- 
tinuum polarization of ~ 15% at a PA of 130°, which is consistent 
with the PA measured by Haas et al. (1992) in optical continuum 
polarization measurements, and also with with the large scale dusty 
disk seen using imaging polarimetry (Minchin et al. 1991). In the 
polarisation triplot, there appears evidence of changes in line polar- 
isation, in particular a rotation in the PA. This is most likely a real 
feature because the emission line is resolved in our data while the 
feature is larger than the formal errorbars in both the PA spectrum, 
and the QU graph. 

We note that this effect may well be similar to the line ef- 
fects measured in many Herbig Ae stars by Vink et al. (2002). In 
their work on these objects they found that the polarisation changes 
across Ha were narrower than the width of the intensity profile it- 
self - which is inconsistent with depolarisation. Instead, the data 
showed line polarisations that are thought to be due to Ha emission 



originating from a compact source of photons that scatters off an 
exterior rotating medium. In contrast to the effect of depolarisation 
observed for MWC 349A, these line polarisations may yield kine- 
matic information. The PA rotations across the line profile translate 
into "loops" when the same data are plotted in QU space. 

Indeed, the excursion in QU space resembles a loop, rather 
than the straight line observed for MWC 349A. Comparison with 
models by Vink et al (2005a, see also Wood, Brown & Fox 1993) 
indicates that such behaviour is the result of the presence of a rotat- 
ing disk with an inner hole. For comparison, rotating disks without 
an inner hole, extending down to the stellar photosphere would dis- 
play a double loop in the QU diagram. If the detection is can be 
confirmed, it provides strong clues that GL 490 is surrounded by a 
disk and is still actively accreting material. 



4 SUMMARY 

We have presented the first hydrogen recombination line spectropo- 
larimetric observations in the near-infrared. The results from this 
study are encouraging. Our targets were two optically bright Her- 
big Be type objects (MWC 349A and MWC 342), an optically 
faint massive Young Stellar Object, GL 490 and for comparison 
the well-known classical Be star £ Tau, known to be surrounded by 
an ionised disk. 

Three of these stars show a line effect. For £ Tau and MWC 
349A the line-effect is explained as being due to optically-thin elec- 
tron scattering in a circumstellar disk. These objects have high- 
resolution imaging data available and in both cases the position 
angle of the (larger scale) disk on the sky is consistent with the 
angle derived from the polarization data. This validates spectropo- 
larimetry as a good means to detect disks on small scales. 

The preliminary detection of a PA rotation in GL 490 may 
indicate the presence of a small scale rotating accretion disk with 
an inner hole similar to those recently discovered in Ha in Herbig 
Ae and T Tauri stars (Vink et al. 2002, 2005b). 
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Figure 4. As the previous figures. The polarization data of GL 490, rebinned to 0.25% per pixel. 
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